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Abstract. Recent results of the Fayans energy density functional (EDF) for spherical nuclei are reviewed. A 
comparison is made with predictions of several Skyrme EDFs. The charge radii and characteristics of the first 
2 * excitations in semi-magic nuclei are briefly discussed. The single-particle spectra of doubly magic nuclei 
are considered in more detail. The phonon-particle coupling effects are analyzed including the tadpole term. 


1 Introduction 

In this century the Hartree-Fock (HF) method with 
Skyrme forces yj] dominates in the theoretical description 
of ground states of heavy nuclei. This method is often 
considered as a version of the Energy Density Functional 
(EDF) method by Kohn and Sham yfl, which is based on 
the Hohenberg-Kohn theorem l^l, stating that the ground 
state energy Eq of a Fermi system ia a functional of its 
density p. The Skyrme HF ansatz leads to a rather sim¬ 
ple form of the nuclear EDF. The same is true for the HF 
method with Gogny force |0l which becomes again pop¬ 
ular last decade. For completeness, we mention also the 
relativistic mean field (RMF) model, see recent Ref. ||3l 
and links there, and a newly developed self-consistent ap¬ 
proach, known as a method BCPM (Barcelona - Catania - 
Paris - Madrid) 

We use an alternative EDF method developed by 
Fayans with coauthors Iffl] with more sophisticated 
density dependence. The main, in-volume term of the 
Fayans EDF can be symbolically written as 

= (1) 

2 1 - 1 - 

where p(r) is the total nuclear density, a,a,cr,y being pa¬ 
rameters. The corresponding term of the Skyrme EDF cor¬ 
responds to y-O in this relation. The use of the bare mass 
m*-m is another peculiarity of the Fayans EDF. Both the 
features of the Fayans method are closely related to the 
self-consistent Theory of Finite Fermi Systems (TFFS) 
and represent, in a hidden form, energy dependence ef¬ 
fects inherent to this approach. The latter is based on the 
general principles of the TFFS il2|] w ith inclusion of the 
TFFS self-consistency relations 11 1311 . The final version 
of this approach llldl 1 151] was formulated in terms of the 
quasiparticle Lagrangian which is constructed to pro¬ 
duce the quasiparticle mass operator 2^(r, e). By def¬ 


inition, the latter coincides with the exact mass operator 
£(r, s) at the Fermi surface. In the mixed coordinate- 
momentum representation it depends linearl y on the mo¬ 
mentum squared k^ and the energy s as well il2|] . 

In the TFFS, the effective mass is a product -m*^ 

of the “k-mass” and the “E-mass”. The two effects com¬ 
pensate each other almost exactly Q resulting in m* ^ 1 
which justifies the Eayans choice of the bare mass instead 
of the effective one. The EDF of the self-consistent TFFS 
is found from the Lagrangian X.g according the canonical 
rules. It includes implicitly the Z-factor 


Z(r) = 1 


de 


( 2 ) 
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where the index 0 means that the energy and momentum 
variables are taken at the Fermi surface. Its density depen¬ 
dence can be found explicitly OlSll : 

Z(r) = 2(l+ ^l-4CoTo2p(r)/4) ‘, (3) 

where Co = (dnlds'^y' = is the usual TFFS nor¬ 

malization factor, inverse density of states at the Fermi sur¬ 
face, and the dimensionless parameter Aq 2 determines the 
Z-factor of nuclear matter Zq. In the result, the density 
dependence of the EDF becomes rather complicated Q. 
Fayans with coauthors found iQ] that, in a wide density re¬ 
gion, it can be approximated with high accuracy by a more 
simple expression (ID- Thus, the Fayans EDF method can 
be interpreted as a version of the self-consistent TFFS. 

Till this year. Ref. iflhl] , all self-consistent calculations 
with Fayans functionals were carried out for spherical nu¬ 
clei only. It proved out successful in systematic descrip¬ 
tion of nuclear magnetic ifr^fT^ and quadrupole | 2T ] 
moments and nuclear radii 11221] as well. In the latter case, 
agreement with the data is better than that in all SHF cal¬ 
culations we know. There are two examples more where 
we may compare the two approaches directly. The ener¬ 
gies and B{E2) values for the first 2^ excitations in semi- 
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Figure 1. Charge radii in lead isotopes. Solid line shows the 
Liquid Drop model predictions. 



Figure 2. Excitation energies n)(2|) for tin is otop es. Predictions 
for the SkM* and SLy4 force are taken from 1^ . 


magic nuclei E 13 is one of them. They are described 
with the Fayans EDF much better than in the known SHF 
calculations with SkM* and SLy4 EDFs |24|. The single¬ 
particle energies (SPEs) in seven doubly magic nuclei, 
for which the experimental spectra are known ll25ll . is an¬ 
other example ll26ll . The PC corrections to were found 
self-consistently with account for the tadpole diagram. A 
systematic comparison to the SHE predictions with a pop¬ 
ular HEB-17 EDF 13 was carried out. Even at the level 
of the mean field theory, the Eayans EDE results are signif¬ 
icantly better. Inclusion of the PC corrections makes the 
agreement better yet. 

The article is organized as follows. Sect. 2 contains 
a brief comparison with several Skyrme EDEs of predic¬ 
tions for the charge radii and characteristics of the first 2* 
excitations in semi-magic nuclei. Sect. 3 is devoted to 
the description of SPEs in doubly magic nuclei. Sect. 4 
contains conclusions. 


2 Charge radii and characteristics of 2 | 
levels in semi-magic nuclei 


In Ref. Il22h . systematic self-consistent calculations of 
charge radii /?ch were made on the base of the DE3-a EDE. 
Deformed nuclei were also included into analysis, with 
an approximate taking into account the deformation ef¬ 
fect. Agreement with the data on the level of 0.01 fm 
was achieved, noticeably better than for SLy4 and HEB- 
17 Skyrme EDEs taken for comparison. The HEB-17 pre¬ 
dictions were taken from 12811. whereas for the SLy4 EDE 
calculations were made in 1221] . As an example, the lead 
charge radii are displayed in Eig. 1. We see that the Eayans 
EDE, indeed, describes the radii perfectly well. The HEB- 
17 one reproduces the data reasonably for heavy Pb iso¬ 
topes but fails systematically for A<190, the disagreement 
reaching 0.1 fm. The reason for this is that this EDE er¬ 
roneously predicts rather strong stable deformation for the 
light Pb isotopes which leads to a significant increase of 
Rch values. This problem is discussed in detail in ifn . The 
SLy4 Rch values are systematically higher than the exper¬ 
imental ones at approximately 0.03 fm. This is a typical 
scale of accuracy in describing nuclear radii for different 
Skyrme EDEs. 


A high accuracy of reproducing the charge radii may 
be linked to more adequate density dependence of the 
Eayans EDE compared to the Skyrme one. Indeed, if we 
denote the average error in describing the binding energies 
as 6E and that for the charge radii as dRch, these quantities 
should be, due to the Hohenberg-Kohn theorem ill], pro¬ 
portional to each other, 

(4) 

where the coefficient a depends on the functional we use. 
This almost obvious relation can be proved with such sim¬ 
ple consideration. Let we have an exact EDE £o(p) and 
add to it a small addendum 6&q - Af(p), where the co¬ 
efficient is small, T « 1. Obviously, we get 6E = 
Ajdrf(p(r)) = Aa; then the change of the mean field 
dU - 6fl6p is proportional to A. Let us find the change 
of the density p in the first order of the perturbation the¬ 
ory 5U. It can be easily seen that 6p and the correspond¬ 
ing change of the radius are also proportional to A, 
SRcb — Ab. In the result, we obtain 0 with a - bja. 
As a rule, a fine tuning of the EDE parameters is per¬ 
forming by focusing mainly on reproduction of the nuclear 
masses within a minimal value of 6E. In this case, the ac¬ 
curacy of reproducing the charge radii is proportional to 
the coefficient a. As the analysis of 1 I 22 I 1 showed, for the 
Eayans EDE this coefficient is less than those of the HEB- 
17 and SLy4 functionals. This observation may be linked 
to more sophisticated density dependence of Eayans func¬ 
tional, which allows to incorporate implicitly the energy 
dependence effects. 

In Ref. d, excitation energies 0 J 2 and transition 
probabilities B{E2) of the first 2^ excitations in even tin 
and lead isotopes as well as the quadrupole moments of 
odd neighbors of these isotopes were calculated within the 
self-consistent TEES with the use of the DE3-a EDE. The 
effect of the density dependence of the effective pairing 
interaction was analyzed by comparing results obtained 
by solving the QRPA equations with volume and surface 
pairing. The effect was found to be noticeable giving ev¬ 
idences in favor of the surface pairing. Eor example, the 
2^-energies are systematically higher at 200-400 keV for 
the volume paring as compared with the surface pairing 
case, the latter being on average better. Obtained results 
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Figure 3. B(E2, up) values for tin isotopes. Predictions for the 
SkM* and SLy4 force are taken from 


were compared with predictions of Il24l] for the Skyrme 
EDFs SkM* and SLy4. For tin isotopes, the comparison 
of a >2 and B{E2) values is given in Figs. 2 and 3, cor¬ 
respondingly. On the average, predictions of the Fayans 
EOF with both models for pairing reasonably agrees with 
the data perfectly well, significantly better than those for 
both the Skyrme EDFs. Among the latter, the SkM* EDF 
turns out to be more successful. For SFy4 EDF, the ' '^Sn 
nucleus and its neighbors are very close to the point of the 
quadrupole instability, in contradiction with the data. 


3 Single-particle spectra of magic nuclei 


The major part of modern calculations with the Fayans 
functional in spherical nuclei are made with the EDF DF3- 
a. This is a version 129] of the EDF DF3 SQ with modi¬ 
fied spin-dependent parameters, the spin-orbit k, k' and the 
first spin harmonics gi,g\ which play the role of the effec¬ 
tive tensor forces. The DF3 and DF3-a EDFs both con¬ 
tain 3 non-zero spin-dependent parameters, specifically 
the value of g\-Q is putted. The bulk of the data 1251] for 
SPEs contains 65 spin-orbit differences which permits to 
try to find an optimal set of the spin-dependent parameters 
of the EDF. Such a set DF3-b, with g\+Q, was found in 
but the overall description of these differences turned 
out only a bit better than for DF3 or DF3-a sets. Moreover, 
the average accuracy in reproducing SPEs see Table 1, 
for DF3-b EDF is the same as for DF3 and only a little 
better than for the DF3-a EDF. The latter was chosen for 
systematic calculations as far as it previously proved to be 
successful in description of different nuclear phenomena 
i T^ tBl 2^. For comparison, we found the SPEs of all nu¬ 
clei under consideration for the Skyrme EDF HFB-17 
which is a record-holder in reproducing nuclear masses. 

Figs. 4 and 5 show the comparison to the data 1251] of 
predictions of all the EDFs under discussion for the SPEs 
in ^°®Pb for neutrons and protons, correspondingly. We 
see that both HFB-17 spectra, especially the neutron one, 
are too expanded. This is the result of influence of the ef¬ 
fective mass TO* < TO. The analysis in Ref. S confirmed 
a preference of the choice of the bare mass for describing 
SPEs with Skyrme EDFs without PC corrections. 


Table 1. Average deviations (fc^)rms (MeV) of the theory 
predictions for the single-particle energies from the 
experimental values for magic nuclei. 


Nucleus 

N 

DF3-b 

DF3-a 

DF3 

HFB17 

“^'’Ca 

14 

1.08 

1.25 

1.35 

1.64 

48Ca 

12 

0.89 

1.00 

1.01 

1.70 

56Ni 

14 

1.00 

0.97 

0.85 

1.40 

7SNi 

11 

1.24 

1.41 

1.09 

1.32 

'“Sn 

13 

1.09 

1.17 

1.01 

1.56 

'^^Sn 

17 

0.58 

0.66 

0.55 

1.15 

208 

24 

0.44 

0.51 

0.43 

1.15 

Total 

105 

0.89 

0.98 

0.89 

1.40 


“°Pb Exp. DF3-b DF3-a DF3 HFB17 



Figure 4. Neutron single-particle levels in Experimental 
data from 



Figure 5. Proton single-particle levels in ^°*Pb. Experimental 
data from fj^l . 



Figure 6. PC corrections to the mass operator. The gray circle 
denotes the “tadpole” term. 









































EPJ Web of Conferences 



2 4 6 8 10 

r(fm) 


Figure 7. The vertex gi for the 3, state in 


Let us go to the PC contributions to SPEs. Accounting 
for PC effects, the equation for SPEs and cotTesponding 
wave functions can be written as 




(5) 


where Hq is the quasiparticle Hamiltonian with the spec¬ 
trum e® and is the PC correction to the quasiparticle 
mass operator. After expanding this term in the vicinity of 
e - one hnds 

(6) 

with obvious notation. Here denotes the Z-factor due 
to the PC effects, i.e. that found from Eq. © with substi¬ 
tution of instead of the main mass operator E(e). 

Remember that in the TEES the corresponding Z-factor is 
included in the quasiparticle Hamiltonian Hq. 

The PC correction to the mass operator is displayed in 
Pig. 6 , where gi is the vertex for creating the L-phonon. 
In magic nuclei, it obeys the equation El 


gdoj) = TA(w)gz.(w), (7) 

where A(a;) = f G(s + ujT) G(e - ujT) delilm) is 

the particle-hole propagator, G(e) being the one-particle 
Green function. In obvious symbolic notation, the pole di¬ 
agram corresponds to = (g^, DiGgi), where Di{co) 

is the phonon D-function. 

All the low-lying phonons we consider have natural 
parity. In this case, the vertex gi possesses even T -parity. 
It is a sum of two components with spins S - 0 and 5 = 1, 
respectively. 


gL = gLoir)TLio(n,a) + gLi(r)TLLi(n,a), ( 8 ) 

where Tjis stand for the usual spin-angular tensor opera¬ 
tors. The operators Tuq and Tm have opposite T -parities, 
hence the spin component should be the odd function of 
the excitation energy, gi\ cx cj^. Por the ghost dipole, 
L - I and coi - 0, Eq. O, due to the TEES self- 
consistency relation El . has the exact solution 

0 i(r) = ai(£/t/(r)/£/r)TiM(n), (9) 

where ai - 1/ V2tuZ?i , Bi = HnAjAn is the Bohr- 
Mottelson (BM) mass coefficient El and U(r) is the cen¬ 
tral part of the mean-held potential generated by the en¬ 
ergy functional. 


The second, tadpole, term in Pig. 6 is 


= / ^i^LgiDdcj), ( 10 ) 

where digi can be found El by variation of Eq. 0 in 
the field of the L-phonon: 

Sigi = diTAicjdgi + TdiAicjdgL 

+ TAiuddLgi. (11) 


The quantity 6lA can be readily obtained by variation 
of each Green function in the particle-hole propagator A in 
held gi induced by the L-phonon. The explicit expression 
for the variation di'T can not be found within the TEES as 
in this approach the Landau-Migdal amplitude T is intro¬ 
duced as a phenomenological quantity. In El the ansatz 
was proposed. 




sr(p) 

6p 


6pL, 


( 12 ) 


where 

6pL=Aigi (13) 


is the transition density for excitation of the L-phonon. 
The phonon D-function appears in Eq. (fTOl i after connect¬ 
ing two wavy phonon ends in Eq. (HUl- This corresponds 
to averaging of the product of two boson (phonon) oper¬ 
ators B^Bi over the ground state of the nucleus with no 
phonons. 

The L-phonon excitation energies col and creation am¬ 
plitudes gid) were found by solving Eq. © in a self- 
consistent way with the DP3-a functional. In more de¬ 
tail, the procedure is described in 01911 . All the L-phonons 
we consider are the surface vibrations which belong to the 
Goldstone mode corresponding to the spontaneous break¬ 
ing of the translation symmetry in nuclei El- The coordi¬ 
nate form of their creation amplitudes gdr) is very close 
to that, Eq. 0 , for the ghost 1 phonon, which is the 
lowest energy member of this mode; 


gdr) = addUjdr) +XL{r), (14) 


where the in-volume corTection;\fi(r) is rather small. 

The smallness of the in-volume component xl is 
demonstrated in Eig. 7 for the 3j^ state in which is 

the most collective one among the surface vibrations and 
plays the main role in PC corrections for this nucleus. The 
small spin components 5 = 1 are also displayed. To make 
them distinguishable, they are multiplied by the factor of 
10. The smallness of the spin components is typical for 
L-phonons with a high collectivity. The hrst, surface term 
on the right-hand sight of Eq. (fTO c orresponds to the BM 
model for the surface vibrations Bill , the amplitude ai be¬ 
ing related to the dimensionless BM amplitude Pi as fol¬ 
lows: ul = RPl, where R - roA'^^ is the nucleus radius, 
and rQ - 1.2 fm. 

If one neglects in-volume contributions, the tadpole 
PC term (fTOl l can be reduced to a very simple form; 




tad 

L 


2L + 1 


At/(r). 


2 3 


(15) 
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Table 2. Pole and tadpole contributions (MeV) to PC 
corrections from 3“-states to neutron SPEs in ^“^Pb. 


A 



Ssa 

3d3/2 

-0.150 

0.012 

-0.137 

2^7/2 

-0.142 

0.061 

-0.081 

4si/2 

-0.134 

0.016 

-0.118 

3ds/2 

-0.147 

0.023 

-0.124 

I715/2 

-0.708 

0.204 

-0.504 

1*11/2 

-0.058 

0.198 

0.140 

2g'9/2 

-0.244 

0.076 

-0.167 

3pi/2 

-0.220 

0.053 

-0.167 

2/5/2 

-0.186 

0.094 

-0.092 

^P3I2 

-0.205 

0.056 

-0.149 

1*13/2 

0.057 

0.211 

0.269 

2/7/2 

0.724 

0.091 

0.815 

lhg /2 

-0.014 

0.197 

0.184 


Table 3. Pole and tadpole contributions (MeV) to PC 
corrections from 3“-states to proton SPEs in ^®*Pb. 


A 



6sa 

3pi/2 

-0.375 

0.153 

-0.222 

3p3/2 

-0.371 

0.152 

-0.219 

2/5/2 

-0.278 

0.168 

-0.110 

1*13/2 

-0.534 

0.266 

-0.268 

'^fin 

-0.409 

0.168 

-0.240 

lhg /2 

-0.054 

0.222 

0.168 

3*1/2 

-0.310 

0.143 

-0.167 

2*f3/2 

-0.241 

0.146 

-0.095 

l/*ll/2 

-0.017 

0.246 

0.229 

'2-ds/2 

0.435 

0.147 

0.582 

107/2 

-0.271 

0.197 

-0.074 


The general consideration of the tadpole term within the 
not self-consistent TFFS, with solving Eq. (fTTT i. was car¬ 
ried out in It was found that the in-volume correc¬ 
tions to Eq. (O are, indeed, small for heavy nuclei, e.g., 
for At the same time, for light nuclei, e.g., 

the accurate solution of Eq. (HB diminishes the approxi¬ 
mate value (fTSl) for the tadpole term by ^ 30%. 

Eollowing to lE3] . we neglect the in-volume correc¬ 
tions for all nuclei considered. To find the phonon ampli¬ 
tudes ai, we use the definition - gY^^^KdUldrY’'^™, 
with obvious notation. It should be noted that the values 
of and cr^ are always very close to each other and to 
that which follows from the BM model formula for B(EL): 
B(EL)bm - (3Z/47r)^j8^f?^^ llil]] . where the dimensionless 
BM phonon creation amplitude jdi is related to that used by 
us as ai - Pi^Rj V2L -i- 1, R - 1.2 Eor example, for 
the 3j^ state in we have: - 0.32 fm, - 0.33 fm, 
and - 0.30 fm. 

A comparison of the pole and tadpole PC corrections 
to neutron and proton SPEs, induced by the 3^ state in 

are given in tables 2 and 3, correspondingly. The tad¬ 
pole term is always positive, whereas the pole one is, as a 
rule, negative, and the two contributions are of the oppo¬ 


Table 4. PC effect on average deviations ((5ea)rms (MeV) of the 
theory predictions for SPEs from the experimental values for the 
DE3-a functional. 


Nucleus 

N 

DE3-aH-PC 

DE3-a 

'‘°Ca 

14 

1.30 

1.25 

4SCa 

12 

1.05 

1.00 

56Ni 

14 

0.98 

0.97 

78Ni 

11 

1.34 

1.41 

'0«Sn 

13 

1.21 

1.17 

i32Sn 

17 

0.63 

0.66 

208 

24 

0.38 

0.51 

total 

105 

0.97 

0.98 


site sign. The magnitude of the tadpole term is, as a rule, 
less than the pole one, but comparable with the latter. Es¬ 
pecially it this true for protons. Therefore the sum is often 
essentially less than the pole term alone. A typical sup¬ 
pression of the pole contribution is of 30-50%, but there 
are cases of a stronger suppression, e.g. 2di/2 and lg7/2 
proton states. Moreover, there are several cases than the 
tadpole term dominates: the neutron liii /2 and lhg /2 states 
and the proton lhg /2 and l/iii /2 ones. In these cases, the 
total correction is of the opposite sign as compared with 
the pole term. The general conclusion is that one overes¬ 
timates the PC correction to SPEs neglecting the tadpole 
term, and often it is better to omit it completely than con¬ 
sider the pole term alone. 

The total PC effect to SPEs of all doubly magic nu¬ 
clei under consideration is shown in table 4. We see that it 
makes the agreement essentially better for ^ 

less good in the light nuclei. The latter is explained mainly 
with the use the approximation (fTSl l for the tadpole term. 
According to the exact consideration for "^°’"^^Ca di¬ 
minishes the tadpole term by approximately 30%. Such 
a correction should lead to a better agreement with the 
experimental SPEs. Let us stress once more a high ac¬ 
curacy of reproducing the experimental SPEs in ^°®Pb: 
(fe.i)rms-0-38 (MeV). Indeed, the corresponding result 
1 331] . obtained within the RME approach, is (&2)rms=0.85 
MeV. 


4 Conclusions 

Recent studies with the Eayans EDE for spherical nu¬ 
clei are reviewed and compared to predictions of several 
Skyrme EDEs. The charge radii and the characteristics of 
the first 2^ excitations in semi-magic nuclei are briefly dis¬ 
cussed. The accuracy of reproducing the charge radii is of 
the order of 0.01 fm, which is significantly better than that 
obtained with HEB-17 and SLy4 EDEs. The excitation en¬ 
ergies and B(E2) values of the first 2^ states in the lead 
and tin isotopic chains also agree with the data much bet¬ 
ter than the results |24] obtained with the SLy4 and SKM* 
EDEs. 

The SPEs of doubly magic nuclei are considered in 
more detail. The PC effects in SPEs are analyzed particu- 
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larly. In addition to the usual pole diagram, the consider¬ 
ation includes the so-called tadpole term which is usually 
ignored 113311341] . The latter is considered approximately, 
with neglecting the in-volume components of the vertex 
giir) of creating a surface L-phonon. This approximation 
works well in heavy nuclei, leading to a simple formula 
for the tadpole contribution which can be easily included 
into the calculation scheme. The tadpole contribution is 
often comparable with that of the pole diagram. As a rule, 
these two contributions to SPEs have different signs, their 
sum being often significantly less than that of the pole 
term alone. For the nucleus, the calculation with¬ 
out PC corrections results in the average deviation ((5ei)rms 
from the experimental values equal to 0.51 MeV for the 
Fayans EOF, whereas it is 1.15 MeV for the HFB-17 EDF 
With the PC corrections, including the tadpole contribu¬ 
tion, we obtained very high accuracy: (dei)nns=0.34 MeV. 
For a comparison, the corresponding quantity, found in 
within the RMF approach with PC corrections with¬ 
out the tadpole term, is {(5e.i)rnis=0.85 MeV. 

To conclude, the Fayans EDF leads to a better agree¬ 
ment with the experimental data in all phenomena in 
spherical nuclei considered than several popular Skyrme 
EDFs chosen for a comparison. We relate it to the pe¬ 
culiarities of the Fayans EDF which reflect, in a hidden 
form, the energy dependence effects inherent to the self- 
consistent TFFS. 
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